Ray-tracing in four and higher dimensional black hole spacetimes: 

An analytical approximation 
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We study null rays propagation in a spacetime of static Schwarzschild-Tangherlini black holes in 
arbitrary number of dimensions. We focus on the bending angle and the retarded time delay for 
rays emitted in the vicinity of a black hole and propagating to the infinity. We obtain an analytic 
expression in terms of elementary functions which approximate the bending angle and time delay 
in these spacetimes with high accuracy. We analyze the relative error of the developed analytic 
approximations and show that it is quite small in the complete domain of the parameter space 
for the rays reaching the infinity and for different number of the spacetime dimensions. Possible 
applications of the obtained results are briefly discussed. 
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I. INTRODUCTION 

Effects predicted by the general relativity are impor- 
tant in the vicinity of compact objects, such as neutron 
stars and black holes. One can study such effects by ob- 
serving electromagnetic or gravitational radiation from 
these objects. The first light detected from regions close 
to the black holes was discovered by the ASCA satel- 
lite, [HQ- Now there is much evidence that interesting 
astrophysical effects are connected with the physics in 
the vicinity of black holes. Broadening of X-ray emission 
lines from accretion disks, 0, 0, IE IE 0, H, M, El , and X- 
ray flares, [HI, [12], or quasi-periodic oscillations, (l3lll4j. 
are examples of these effects. To explain these obser- 
vations one needs not only to develop models of these 
phenomena but also to solve equations for the light prop- 
agation in the gravitational field of these objects. 

If the wave-length of the radiation is much smaller 
than the characteristic scale [r g , the gravitational radius) 
one can use the geometric optics approximation and re- 
duce the problem to the study of null ray propagation in 
curved spacetime. Two related problems are of interest: 
(1) how does a distant observer see objects located far 
'behind' a black hole distorting images; (2) how is the 
light emitted by matter in the region of a strong grav- 
itational field seen by a distant observer. For example, 
the radiating matter can be the surface of a collapsing 
body, or the surface of a neutron star, or an accretion 
disk and so on. A similar problem is of potential inter- 
est for higher dimensional black holes. In brane models 
with large extra dimensions mini black holes can play a 
role of 'probes' of extra dimensions. Scattering or emis- 
sion of light or gravitons by such black holes might be 
interesting in this connection (a recent review on higher 
dimensional black holes can be found in [IH, see also 
references therein). 
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There exist a lot of publications where ray-tracing in 
Schwarzschild and Kerr geometries has been discussed in 
detail. It is well known that the Hamilton-Jacobi equa- 
tions for light rays in these geometries allow separation 
of variables. This property is connected with the exis- 
tence of an additional quadratic in momentum integral 
of motion associated with the Killing tensor, pjj] ■ In 
the Schwarzschild 4D geometry expressions for scatter- 
ing data (such as a bending angle and time dela y) c an be 
written explicitly in terms of elliptic integrals, |17j . [l8j |. 
In the Kerr metric such quantities can be expressed in 
terms of generalized hyper-geometric functions, jl9j |. 

Recently it was shown that higher dimensional rotat- 
ing black holes in many aspects are similar to the 4D 
Kerr black holes. Namely, the most general Kerr-NUT- 
(A)dS metric always possesses a so-called principal con- 
formal Killing- Yano tensor [20l . [2lj whi ch g enerates a set 
of the second rank Killing tensors [HI, [23|. As a result 
the geodesic equations in such spaces are completely inte- 
grable [IS 01, an d the Hamilton-Jacobi, Klein-Gordon 
and Dirac equations allow complete separation of vari- 
ables [HI, [H[. It means, that a solution of the geodesic 
equations can be obtained in quadratures. However, even 
for the non-rotating black holes in higher dimensions a 
solution cannot be expressed in terms of known special 
functions (for a general discussion of hidden symmetries 
and separation of variables in higher dimensional black 
holes see recent reviews [13, [Hj] and references therein) . 

For scattering at small angles the integrals arising in 
these problems can be estimated by using the perturba- 
tion theory (see e.g. [11]). One can use also numerical 
calculations, as it was done in [30| where the capture 
cross-sections for a five-dimensional rotating black hole 
was studied. Usually the ray-tracing, used for the recon- 
struction of spectrum and light curves as seen at infinity 
for light emitting patterns, requires repeating calcula- 
tions for very many rays. To reduce the time and cost 
of the calculations, and to be able to study analytically 
different qualitative observed features of the radiation it 
is useful to have simple analytical expressions, for exam- 
ple in terms of elementary functions, which approximate 
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accurately enough the scattering data. In this paper we 
develop a scheme for obtaining such approximations for 
four and higher dimensional non-rotating black holes. 

An elegant analytical formula approximating bending 
angles for null rays in the 4D Schwarzschild geometry 
was proposed by Beloborodov, The accuracy of 

the Beloborodov's approximation for light rays passing 
the gravitating object at R = 3r g is of the order of 1%, 
but unfortunately for the rays passing closer to the black 
hole its accuracy becomes worse and reached, e.g., 10% 
for those ra ys which pass at R = 2r g . In the more re- 
cent paper [32j there was proposed a better analytical 
approximation for the bending angle and time delay . It 
has an accuracy of 2 — 3% for rays emitted at the radius 
R > 2r 9 . However, for rays emitted inside R — AM the 
accuracy of this approximation diminishes. The reason is 
this: the bending angle for rays with impact parameter 
close to the critical value A* = 3\/3M passing close to the 
critical radius r* = 3M becomes (logarithmically) large. 
The asymptotic behavior of the bending angle for near 
critical rays was studied in [33[ . Very recently an analyt- 
ical procedure for studying gravitational lensing in the 
strong deflection regime was proposed, (34|. In essence 
[34j extracts from the bending the angle the divergent 
term and studies the behaviour of the deflection of light 
near the critical radius. As a result [3~ij constructs an ap- 
proximation to the four dimensional gravitational lensing 
equation in the strong deflection limit for the Kerr black 
hole and spacetimes which have a stationary spherically 
symmetric line element. In the Schwarzschild case for 
photons emitted at R = AM and scattered to infinity at 
R = 3.05M, for example, the method outlined in [34| 
results in a relative error for the bending angle on the 
order of 0.5% when the impact parameter has the value 
5.198M, and for those photons escaping to infinity from 
inside the critical radius at R = 2.5M, for example, the 
relative error is of order 0.05% when the impact param- 
eter is 5.144M. 

In the present work we propose an improved analyti- 
cal approximation for the ray tracing problem in four and 
higher dimensional static black hole space-times. We an- 
alyze the following problem: suppose a ray is emitted 
at the radius ro with impact parameter A. We obtain 
an analytic expression for the bending angle and time 
delay for such rays which is uniformly valid for the two 
parameter set {ro, A} specifying the ray. To make an ap- 
proximation possible we first extract from the integrals 
for scattering data the contributions which are logarith- 
mically divergent near the critical trajectories. The key 
observation of the present work is that after this proce- 
dure the remaining part can be approximated with very 
high accuracy by a function of one variable. By finding 
a proper approximation for this function one can obtain 
a very accurate uniform approximation for the required 
quantities. 

The paper is organized as follows. In section II we col- 
lect formulas for the ray propagation in the Schwarzschild 
and Tangherlini geometries. In sections III and IV we 



derive the expressions which we use to approximate the 
bending angle and time delay, and study the errors of 
the approximation. Section V contains a discussion of 
the obtained results and their possible applications. 



II. NULL RAYS IN A STATIC BLACK HOLE 
GEOMETRY 

A. Basic equations 

We consider null rays propagating in the background 
of a Z?-dimensional Schwarzschild- Tangherlini metric 



ds 2 = -fdt 2 + f^dr 2 + r 2 dn 2 n+1 , 
where n = D — 3, 

16ttM 



/ = i- - 



(n + l)A n+1 



(1) 



(2) 



Here M is the mass of the black hole and dtt 2 l+1 is a line 
element on a unit (n + l)-dimensional sphere, S n+1 , 

dfi? = de 2 a = d^ 2 , d ^l+i = do l + sin2 , (3) 



6i>i e [0, tt] , = 6 € [0,2tt], 



and A n+ i is the area of S n+1 , 



A 



71 + 1 



2 7r ( n + 2 )/ 2 

r(2±2) 



(4) 



Here T(z) is the Euler Gamma function. We use units in 
which the D-dimensional gravitational coupling constant 
G and the speed of light c are equal to 1. 

It is easy to show that similar to the D = A case a 
photon trajectory in the spacetime (fT]) lies always within 
a plane. Without loss of generality we assume this plane 
to be the equatorial plane, i.e. 9i = tt/2, i = l,..,n. 
We choose an affine parameter ( so that photon's D- 
momentum is p^ — dx^/d( and for our choice of the 
coordinates we have 



P 



(p*,p r ,p*,0,...,0) 



(5) 



The energy E — —p t and the angular momentum L = 
are integrals of the motion. Since a trajectory with L < 
can be obtained from a trajectory L > by a simple 
reflection <fi — > —<fi we assume that L > 0. Using the 
integrals of motion and the relation p^p^ = one can 
write the equations of motion in the form 



r = a r E Z , 
i = E/f, 



Z=y/1- X 2 f/r 2 
■ = L/r 2 . 



(6) 
(7) 



Here the dot over an expression means its derivative with 
respect to the affine parameter £ and a = ±1. For the 
outward moving photon ay = 1, while for the inward 
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FIG. 1: The above figure shows qualitatively the different 
types of motion for light rays in a spherically symmetric black 
hole background. Line 2 shows a ray which has a turning 
point, A. It corresponds to the scattering problem, that is 
when the ray coming from infinity is directly scattered back 
to infinity. Ray 1 corresponds to a critical trajectory when an 
incoming photon spirals around the circle (the critical radius) 
an infinite number of times. Ray 3 describes a motion of a 
photon coming from infinity which is captured by the black 
hole. 



moving one ay = — 1. A change of the sign of oy occurs 
at a turning point r* defined by the relation 



l-\ 2 .U/rl = Q, f t = f(r 



(8) 



By excluding the affine parameter the equations ©-(HI) 
can be written in the form 



d(j> 
dr 



cr r X 



dt 
dr 



(7 r 

Jz 



(9) 



where A = L/E is the photon's impact parameter. 

Consider a photon emitted at point r = r . We call 
such a photon forward (backward) emitted if oy > 
(a r „ < 0) at the point of emission. For a photon prop- 
agating from r to infinity we define a bending angle as 
A<p = (f>(oo) — 4>{ro). The bending angles for the forward 
and backward emitted photons are, respectively, 



A</>+ = $(r ) , $M 



Xdr 



(10) 



A0_ = 



^=2$(r.)-*(r ). (11) 
r z Z 



Similarly, the time of arrival of a forward emitted pho- 
tons as measured by an observer at radius r± is 



(°) _ Ae) 



dr 

Jz 



(12) 



The integral in (|12| is divergent when the upper limit 
tends to infinity. This divergence reflects a simple fact: 
namely, the time required for a light ray to reach infinity 
is infinitely large. For this reason it is more convenient 
to deal with the retarded time 



u = t — f 



dr 

dr=j 



After simple transformations one has 



(13) 



(14) 



r(r ) 



dr 
1 



Z 



- 1 



» 



» 



+ 2(f - u) + 2T(r*) - T(r ) 



(15) 



(16) 



where fo and f* are tortoise-like coordinates, (|13[) . for 
the point of emission and turning point, respectively. 

In what follows we focus on the functions <& and T since 
the bending angle and retarded time for both forward and 
backward emitted photons can be expressed in terms of 
these quantities. 

In the 4 dimensional spacetime, n = 1, the integrals 
(10) and (15) can be written in terms of elliptic func- 
tions. In the higher dimensional case the integrals can- 
not be written in terms of known special functions. Our 
aim is to study these objects in a spacetime with arbi- 
trary number of dimensions as functions of two variables, 
the impact parameter A specifying the photon trajectory, 
and the point of emission, ro- It should be emphasized 
that knowledge of these scattering data allows one to ob- 
tain expressions for other quantities which might be of 
interest. For example, if one wants to know what are the 
bending angle and retarded time delay for a photon with 
given impact parameter propagating from one point, ro, 
to another, r%, it is sufficient to calculate the differences 
between the corresponding quantities A<fr\(ri) — A0^(ro) 
and Au A (ri) - Aux(r ). 



B. The equations of motion in a dimensionless form 



The problem contains only one dimensional parame- 
ter r g , which determines the scale. It is convenient to 
introduce dimensionless quantities 



r 9 

Q=—, 
r 



In these variables 



'7o 



s 5s., i = -. 

r r„ 



f = l-q n , Z=y/1-Pq2f. 



(17) 



(18) 



The motion of the photon is possible only in the region 
where 



< I < L 



(19) 



(Notice that we consider only non-negative impact pa- 
rameters.) On a plot in (q, I) variables this region lies 
below the line l m ax(q)- The allowed regions for differ- 
ent number of dimensions (n > 1) are similar. Figure 2 
shows the allowed domain for n = 1. The minimum of 
the function l max (q) is at q = q« , where 



l/n 



(20) 
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FIG. 2: Function l m ax(q). Motion of photons is allowed in 
the region below this curve. Points q — and q = 1 corre- 
spond to the infinity and the horizon, respectively. Horizontal 
lines describe the photon trajectories. Line 2 corresponds to 
a photon scattering, while line 3 corresponds to the photon 
capture. Line 1 is critical. It corresponds to a photon with 
the critical impact parameter U. Domain / is a region below 
= lmax{q) lying to the left from a vertical line 4. A region to 
the right from line 4 and below line 1 is domain II. In what 
follows we consider only photons emitted from domains I and 
//, which can reach the infinity. 



At this point I — Z* and / = /* where 



L = 



n + 2 f n + 2 



/* = 



n + 2 



(21) 



Before going further we distinguish the possible qual- 
itatively different types of the light ray trajectories (see 
Fig. [2|). Rays with I < U can propagate from infinity 
(q = 0) to the horizon (q = 1) or travel to infinity from 
a point q satisfying q < 1. Rays with I = can either 
approach q — q* coming from the black hole (q > q*) 
or from infinity q = 0). In doing so they spiral around 
q = q* an infinite number of times. Rays with I > I* 
have a turning point. They can either propagate from 
infinity to some pcriastron and return to infinity or from 
the black hole horizon to an apastron and after fall back 
to the black hole. The latter case is not of much interest 
for astrophysical and other applications and we shall not 
consider it. 

Using the above introduced dimensionless quantities 
we can rewrite the expressions for the bending angle and 
time delay in the following form 



90 Idq 



T 



1 



90 dq_ 
<ff\Z{q) 
Z\q) = \-l\\\-q n ) 



- 1 



(22) 

(23) 
(24) 



Let us consider the function Z(q) in the interval from 
q = (infinity) to q = 1 (horizon). It monotonically 
decreases from its value 1 at q — until the minimum 



(25) 



at q = g*, and then monotonically increases until it has 
the value 1 at q = 1. It means that one can expect that 
for the integrals in (|22|) and ([23]) the main contribution 
comes either from the region near the point if it enters 
the integration domain, g» € (0,qo), or from the region 
near the end point go hi the opposite case. We use this 
remark to construct an approximation for these integrals. 



III. APPROXIMATING THE BENDING ANGLE 



A. Type I rays 



1. Leading part 



First we study the rays emitted from the domain I 
{type I rays). Let us introduce a new coordinate y related 
to r and q as follows 



y = l-ro/r, 



3o (1 - v) ■ 



and denote 



pa = qo/q*, p = q/q*=po(^-y)- 



p 



Po 



D 



l/k 



(26) 



(27) 
(28) 



By fixing the parameters P and B one uniquely specifies 
a ray and its emission point. For the rays of type I one 
has < P < 1 and < B < 1. 
Using these notations one gets 



/ = 1 



2P 



(1-2/)", 



(29) 



and the expression for the bending angle for the forward 
emitted photons takes the form 



1 vdy_ 

Q 



= vfoz = ^f -^(i~ y rf. 



One also has 



lv=o = (1-S)(1- 



2P 



) 



(30) 
(31) 

(32) 
(33) 



Q 2 \ y =i = / = l-2P/(n + 2). 

We study now $ as a function of (P, B) in the domain 
I : < P, B < 1. Our goal is to obtain an analytic 
expression uniformly approximating $ in this region. 

For the interval (0,qo < q*) in the domain / the func- 
tion Q has its minimal values at y = which corresponds 
to q = q . We denote by Q 2 the expansion of Q 2 in y up 
to the second order 



2P 



Q 2 = a + 2by + cy 2 

a = (1-S)(1 

n + 2 

b = B(l-P), 

c = -B[l - (n + l)P] 



(34) 
(35) 



In the domain I the parameters a and b are always pos- 
itive, while c > for P > l/(n + 1) and c < for 
P < l/(n + 1). At the point y = 0, Q 2 and its first 
two derivatives coincide with the similar quantities for 
the function Q 2 . In particular, at in the vicinity of the 
point (P = B = 1) one has 



Q ~ <9 ~ Vm/ + • ■ • • 

At the point y = 1 one has 

^ n + 2 n + 2 



BP. 



Let us denote 

and present $ in the form 

$ = $ + A$ . 



(36) 



(37) 



(38) 



(39) 



Both integrals, <& and $ are logarithmically divergent at 
the lower limit, y = 0, for the point P = B = 1 of the 
parameter space. This point corresponds to a trajectory 
with the impact parameter I = I* emitted at q = q*. 
As a consequence of the property (|3"6"|) the divergences of 
the both quantities $ and $ are identical and hence A$ 
remains finite at P = B = 1. The integral (f38|) can be 
easily taken. We shall discuss now its exact value and 
later we shall find a suitable approximation for A$ valid 
in the total domain /. 

The form of the expression for $ depends on the sign 
of c and on the sign of the discriminant D (see figure GU, 



D = ac - b 2 . 
The equation of the line where D — is 
B= [(n + l)P- l](n + 2-2P)/(nP(n- 



(40) 

1-P)). (41) 

This line is shown in Fig. [3] It starts at B — and 
P = l/(n+ 1) and ends at the corner P = B = 1. The 
vertical lines in the figure [3] correspond to the fixed value 
of P (and hence of q and r). The lines which start at the 
point (0, 0) in this figure correspond to the fixed values of 
the impact parameter I and are described by the equation 



B = l 2 qlP 2/n [l - 2P/(n + 2)] . 
Inside the domain / for c < one has D < and 
^cB 



(42) 



4> 



D 



'-[nPy/R- (1 -P)VI] . (43) 



The parameters A and D which enter this relation are 
defined by (137]) and (|40|) . respectively. 

Taking the limit c — > 0, (P ^ 0), in this expression one 
obtains 



-(Va + 26- Va) . 



(44) 
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FIG. 3: The above figure shows the domain / of Figure 2 
in (P, B) coordinates. The line 1 originating from the point 
(0,0) and terminating at the point (1, 1) corresponds to pho- 
tons with the critical value of the impact parameter. Lines 
originating from the point (0, 0) which are above line 1 cor- 
respond to line 2 in figure 2 and lines below it correspond to 
line 3 in figure 2. The line c where P = l/(n + 1), separates 
two regions where c < (to the left) and c > (to the right). 
The discriminant D vanishes on the line D. A region with 
negative D is to the left of this line, and the region with pos- 
itive D is to the right of the line D. The plot is shown for 
n = 2. For other values of n the structure is similar. 



Finally for c > one has 



$ = — In 

\/c 



nBP 



<cA 



(45) 



We shall need later the expression of <t for the special 
value P — 1. In this case c > 0. We denote this quantity 
by For this choice one has go = Imax = I* and v 
takes the value /i = i/Z*. Taking this limit in the formula 
(1431) one obtains 



1 



dn 



fi^n + 2+ y/(n + 2)[i 2 + 1 - fi 2 

For a trajectory close to the critical one, /i ~ 1, 

1 



ln(l-M) 



(46) 



(47) 



2. Approximating A«I> 

Now we focus our attention on the quantity A$ which 
describes the difference between the exact value of the 
bending angle $ and its leading part $ which we calcu- 
lated in the previous subsection. We have plotted A$ as 
a function of po and v in the domain /. After a study 
of these plots we found a remarkable fact. Namely, A$, 
which by its definition is a function of 2 variables, P and 
B, can be approximated accurately enough by a function 
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of one variable, which is a linear combination of po and v. 
We found that this fact is valid not only in 4 dimensions, 
but in higher dimensions as well. Figure[3]illustrates this. 
It gives an example of the 3D plot of A<!> constructed by 
Maple for n = 2 (D = 5) with the orientation option for 
the view angles chosen to be [—47°, 90°]. The 2D surface 
in this projection is a slightly broadened curve. Similar 
orientation angles can be found for other dimensions. Us- 
ing this fact and by fitting the corresponding 2D plots, 
one can find the following expression which approximate 
A$ with a very high accuracy 



A$ = -k 9 {y/T- 
.t$ = \/2[cos( — 



(a;* - l) 2 - l)9(x$ - 1) 
I- ip$)po + sin(- + 



(48) 
(49) 



Here G is the Heaviside step function. The parameters 
fc$ and "0* for different number of spacetime dimensions 
are given in Table [J 



Dimension 


fc$ 




Percentage error 5$ 


D=4,n=l 


0.47 


tt/45 


[-0.77% , 1.78%] 


D=5,n=2 


0.49 


7r/90 


[-1.20% , 1.69%] 


D=6,n=3 


0.52 





[-2.10% , 1.71%] 


D=7,n=4 


0.55 





[-2.89% , 2.58%] 



TABLE I: Approximation parameters and V>*, and range 
of relative errors (5$ for <E> in the domain I for the spacetime 
of dimension D = n + 3. 

This observation implies that one can approximate $ 
by the following expression 



Here $ is given by (|4"3"]) -(|45 p . and A<t is given by 
(|49p . The relative error of this approximation is 

$ - $ a A$ - A$ 



(51) 



Such a relative error calculated in the domain / belongs 
to some error interval. These error intervals for the space- 
times with different number of dimensions are given in 
Table HI The figure 5 illustrates how the relative error 
is distributed in the space of parameters (po; v )- This 
particular plot of the surface 5$, is again shown for n = 2. 
The plots for other values of n are similar. It is evident 
from the figure that our approximation for the bending 
angle works very well everywhere in the domain / in- 
cluding the region near (po = 1> v = 1) which is precisely 
where the integral expression for the bending angle di- 
verges. 



B. Type II rays 

1. Leading part 

For the light ray emitted from the domain 7T the main 
contribution to the bending angle is from an interior 
point of the integration domain. For this reason the 
above procedure for obtaining the analytical approxima- 
tion for the bending angle must be slightly modified. We 
rewrite the expression (|2"2"| for the bending angle in the 
form 



$ sa $° = $ + A$ . 



(50) 



4> 



pa 



dp 



(52) 
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where, as earlier p — q/q*, /U = l/l* and 



Z\p) = l_!^ M V(l-^P n ). (53) 

Since po > 1 an d the function Z(p) has its minimum at 
p = 1, it is convenient to split the integration domain in 
(15^)1 into intervals (0, 1) and (l,po). To approximate the 
first integral 



** = A* 



n + 2 f 1 dp 



(54) 



one can use the expression obtained in the previous sec- 
tion 



(55) 



where <&* is given by (|4"o]l . while A$* is A$, given by 
(|48[). calculated for po = 1- Thus, for our purposes, it 
is sufficient to discuss the approximation of the following 
quantity 



2 r po dp 



(56) 



Expanding the function Z 2 (p) in the powers of (p — 1) 
near its minimum, p — 1 = 0, and keeping the terms up 
to the second order one obtains 

Z 2 = (l-p 2 ) + (n + 2)^ 2 (p-l) 2 . (57) 
As earlier, we substitute Z(p) instead of Z(p) in (|56|) 



n + 2 



dp 
i #(p) 



(58) 



This integral can be taken explicitly and defining z = 
H\/n + 2(po — 1) one gets 



» 1 . z + Jz 2 + 1 - fi 2 
$77 = — f=1h / 



(59) 



We now turn our attention to the relative error in the 
approximation of $77 by $77. The relative error is 



S' 1 ' 



12 



i7 



17 



(60) 



The relative errors in the dimensions considered are given 
in Table II and the figure 6 is a plot of Sfj in five dimen- 
sions, the plots for other dimensions are similar. One can 
improve the accuracy of the approximation for $77 as it 
was done for rays emitted in the domain /. However, 
for most of our purposes it is sufficient to approximate 
$77 by $77. This happens because the relative contri- 
bution to the bending angle by the inner part of the ray 
lying within the domain //, (1 > q > q*), is smaller than 
the contribution from the part of this ray in the domain 
/. For this reason the contribution to the 'inaccuracy' 




FIG. 6: This is a plot for n = 2 of the surface 8 U for rays in 
domain II for po satisfying 1 < po < 1.41, < p. < 0.99 when 
it is rotated to have an orientation [-9°, -107°]. We don't show 
the region for 0.99 < fj. < 1 because the plot descends steeply 
to zero in this region and it is difficult to show the nature of 
the plot of this region on this scale. 



of the inner part to the total 'inaccuracy' is suppressed. 
This can be seen from the Table II which contains both 
relative errors, and Sfj. Thus for the approximation 
of the inner contribution $77 to the bending angle it is 
sufficient to use only the leading term $77. 

Using this information we approximate the bending 
angle for light propagating to infinity from domain // by 



$ a = + A$* + $77 



(61) 



The use of our approximation in the previous section and 
$77 accounts for the reasonably good relative error of our 
approximation of the bending angle for light coming from 
domain //. The relative error of our approximation is 

$ — <f>° 

<5* = ■ (62) 

The Table II contains the relative errors appropriate to 
each value of n considered. The figure 7 is a plot of 5$ 
for n = 2. We can see that our approximation remains 
good in the region where the bending angle diverges. 



IV. APPROXIMATING THE TIME DELAY 
A. Type / rays 

1. Leading part 

The time delay ft ' — t^ e \ (fF2"j) . and the correspond- 
ing delay in the retarded time u, (HP, is defined by the 
function T, (|15p . or its dimcnsionless form r 



T = r 9% 1t > r = 



dy 



(1 - v) 2 f\Q 



- 1 



(63) 
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l~)impn Qinn 


Po 


Pprppnta trp prrnr 


Pprppntaffp prrnr 


D=4,n=l 


pa = 1.49 


[-4.38%, 0.00%] 


[-1.30%, 0.83%] 


D=5,n=2 


po = 1.41 


[-5.69%, 0.00%] 


[-1.80% , 1.15%] 


D=6,n=3 


po = 1.35 


[-6.64%, 0.00%] 


[-2.11% , 1.72%] 


D=7,n=4 


Po = 1.31 


[-7.50%, 0.00%] 


[-2.25% , 2.56%] 



TABLE II: Above are the relative errors associated with the 
quantities Sfj and 5$ for tf>u and $ in the domain II for po 
satisfying 1 < po < q* 1 for the spacetimes of dimension D — 
n + 3. Also included are the maximal values of po considered 
for the endpoint of integration. 




0.75 



FIG. 7: This is a plot for n — 2 of the surface <5$ xi for rays in 
domain // for po satisfying 1 < po < 1.41,0 < p < 0.99 when 
it is rotated to have an orientation [-48°, 62°]. We don't show 
the region for 0.99 < p < 1 because the plot descends steeply 
to zero in this region and it is difficult to show the nature of 
the plot of this region on this scale. 



We now construct an analytic approximation for r in 
the same way as it was done for $ in the previous section. 
It is easy to see that in spite of the factor (1 — y) 2 in the 
denominator of the integrand in (|63[) , the integral is finite 
at y = 1. The only singular point of r in the domain / 
is its logarithmic divergence at y = for the limiting 
values of the parameters po = v = 1. To extract this 
divergence we subtract from r a similar integral, where 
instead of Q we use Q given by (f3"4")l . In order to be 
able to obtain the value of the integral in an explicit 
form in terms of elementary functions, as earlier, we make 
another modification of the integrand. Namely, instead 
of H = (1 — y) 2 f we use its quadratic expansion near 
y = of the form 



H 



fo + 2py + ly 2 



(64) 



fo = 1- 

7 



2P 



n + 2 
l-(n + l)P, P 



/3=-(l-P): 

Po- 



(65) 
(66) 



It is easy to see that fo > and that the discriminant 
of the quadratic in y, is of the form 



/o7 - /? = - 



nP(n +1-P) D — f c 



(67) 



and is non-positive. The parameters a, b, and c, which 
enter the expression Q, are related to the coefficients 
in as follows 



(l-^ 2 )/o, 6 



(68) 



Let us emphasize that the structure of the logarithmic di- 
vergence of the exact integral (|63[) at y = is preserved in 
our approximation since H(y = 0) = fo = [(1 — Z/) 2 /] |y=o- 
We define the quantity f as follows 



n H 



To 



Q 



(69) 



We note that H has a root in the interval [0, 1] . Rewriting 
the integrand in (|6"9"]) as follows 



1 

H 



To 



Q 



l 



TQ + Q 2 



(70) 



one can see that it does not result in a divergence. Defin- 
ing the following quantities: 

M = f c-D, 

At = D + cy r hA± (c + b)VM, (71) 
B± = D + cyTa±bVM, 
where, as earlier, D = ac — b 2 , the explicit form of f is 

'A. B 



2VM 



hi 



A-B A 



(72) 



In the case c = 0, (B ^ 0), f has the following form 



In 



, / a + 2b — fo 



2b 

,,2 



a - fo 



— tanh 1 
b 



T(VaT2b - y/g) 
fo 



a-J a + 2b 



(73) 



For the boundary of the domain I where q = q* (po 
1) the expression for f simplifies and takes the form 



VnT2 
t* = — ln(Z+/Z_). 



(74) 



Z ± = Bn^/n + 2 ± ( n — n B + ny/{n+l)B + 1) 
For the near critical rays (// — 1 w 0) 

n + 2 



v/tT+2 
2n 



In(M-l) 



hi 



2 v / rT+2 + ' 



(75) 
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FIG. 8: This is a plot for n — 2 of the surface At when it is 
rotated to have an orientation [-45°, 90°]. We note that we 
use the orientation [-43°, 90°] in our approximation because 
it better approximates a curve. 



This time is logarithmically divergent at fx — > 1. 

As in the case for $ we can express the quantity r as 

t = f + At , (76) 

the quantity At being finite for rays moving along critical 
trajectories. 



Dimension 




IpT 


Percentage error S T 


D=4,n=l 


0.66 


tt/45 


[-1.68% , 2.71%] 


D=5,n=2 


0.55 


tt/90 


[-2.19% , 2.60%] 


D=6,n=3 


0.52 





[-3.28% , 2.31%] 


D=7,n=4 


0.52 


tt/60 


[-1.73% , 2.66%] 



TABLE III: Approximation parameters k T and tp T , and range 
of relative errors S T for t in the domain II for the spacetime 
of dimension D = n + 3. 




FIG. 9: This is a plot for n = 2 of the surface 5 T when it is 
rotated to have an orientation [29°, 46°] 



2. Approximating At 

As earlier, by studying of the plots for At as a func- 
tion of 2 variables po and v one can observe that they can 
be approximated by a function of one variable, which is 
a linear combination of po and v. For example, figure [5] 
shows the 3D plot of At for n — 2 (D = 5) with the orien- 
tation option chosen to be [—43°, 90°]. The 2D surface in 
this projection is again a slightly broadened curve. This 
allows one to approximate At with a very good accuracy 
by a function 

At - -fc^oKv 7 ! ~ (ir - I) 2 - l)9(a: T - 1) , (77) 

x T = \/2[cos(^ + ip T )po + sin(— + VvM • (78) 

The corresponding values of k T and ip T are given in Ta- 
ble IIII1 Thus we can approximate r as follows 

t « r a = f + At . (79) 

The figure 9 is a plot of the surface 5 T for the domain 
J, < {poj^} < 1- Here we have shown S T for n = 2. 



The plots for other values of n are similar. The relative 
error of this approximation is 

T~r a At- At , n . 

5 T = = . (80) 

T T 

The relative errors S T for a spacetime with different num- 
ber of dimensions are given in Table IIIH It is evident 
from the figure that our approximation for the time de- 
lay works very well in the region near (j>q = l,v = 1) 
which is precisely where the integral expression for the 
retarded time diverges. 



B. Type II rays 

1. Leading part 

In this section we investigate the retarded time for light 
rays that proceed to infinity from the domain iT where po 
satisfies 1 < p < ((n + 2)/2) 1 /". As in the case of type I 
rays we express the retarded time for a light ray emitted 
in domain // in terms of dimensionless quantities. In this 
instance we use the quantities po and fi as defined earlier. 
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After reformulating the expression for the retarded time, 
(fT5)) . in terms of these variables we have 



T 



J(pi,Po) 



1 + - J(0,po) 
n , 

dp 



pi Z(p)(l + Z(p)) 



(81) 
(82) 



Since we are in a region where po > 1 arid Z(p) has a 
minimum at p = 1 we again split the integration domain 
in (|8Tj) into two intervals, (0, 1) and (l,p )- In the first 
interval we can write the contribution to the retarded 
time as 



t* = jit 



1 



J(0,1) 



This integral can be approximated by 

Tat "^"^ 7";); _ (~ ^^T^fe . 



(83) 



(84) 



where f* is given by (f?4")) and Af* is Af evaluated at 
Po = I. Therefore in approximating the retarded time in 
domain II it is sufficient for our purposes to consider the 
following integral 



Til = A* 



1 + - J(l )Po ). 

n , 



(85) 



Identically to what we did earlier we substitute Z{p) 
given by (f57|) for Z(p) to obtain 



fix =11 z \ 1 + - ) J(1,po), 



</(l,Po) 



Po 



Pl Z(p)(l + Z(p)) 



(86) 
(87) 



This integral is calculated explicitly using Maple and is 
given by 



v 71 - M 2 + + 1 - A* 2 

tii = — In ■ 



2n " y/l - M 2 + z 2 + i _ _ ^ 
We now focus our attention on the relative error 5 



(88) 



ii 



S T n = 



tii - tii 
th 



th 



(89) 



of the approximation of tjj by f//. In Table IV we 
present the relative errors of our first order approxima- 
tion of th for dimensions 4, 5, 6, 7. The figure 10 is a plot 
of 5jj in five dimensions, the plots in other dimensions 
are similar. As earlier, examining the results presented 
in Table IV we can conclude that tu provides a good 
first order approximation of the function th in the do- 
main //. Therefore in approximating the retarded time 
for light going to infinity from domain II we use 



Af* 



TU ■ 



(90) 



-0.06 G: 




FIG. 10: This is a plot for n = 2 of the surface S\j in domain 
II for po satisfying 1 < po < 1.41, < p, < 0.99 when it is 
rotated to have an orientation [-6°, -106°]. We don't show 
the region for 0.99 < /i < 1 because the plot descends steeply 
to zero in this region and it is difficult to show the nature of 
the plot of this region on this scale. 



and the corresponding relative error of our approximation 
for the retarded time delay is 



(91) 



This relative error for each dimension considered is shown 
in Table IV. The figure 11 is a plot of 8 T for n — 2. We 
can see that our approximation is very good in the region 
fi ~ 1, which is where the retarded time diverges. 



Dimension 


po 


Percentage error 8 T n 


Percentage error S T 


D=4,n=l 


po = 1-49 


[-5.18%, 0.00%] 


[-2.00%, 0.72%] 


D=5,n=2 


Po = 1-41 


[-6.73%, 0.00%] 


[-2.52%, 0.16%] 


D=6,n=3 


po = 1.35 


[-7.83%, 0.00%] 


[-3.24%, 0.67%] 


D=7,n=4 


po = 1.31 


[-8.83%, 0.00%] 


[-2.59%, 2.67%] 



TABLE IV: This table shows the relative errors associated 
with the quantities <5J/ and 5 T for our approximation of tii 
and r in the domain II for po satisfying 1 < po < q* 1 for the 
spacetimes of dimension D — n + 3. Also included are the 
maximal values of po considered for the endpoint of integra- 
tion. 



V. SUMMARY AND DISCUSSION 

In the present paper we have investigated ray-tracing 
problem for photons propagating in four and higher di- 
mensional spherically symmetric black hole backgrounds. 
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FIG. 11: This is a plot for n = 2 of the surface 5 T for rays in 
domain II for po satisfying 1 < po < 1.41,0 < fi < 0.99 when 
it is rotated to have an orientation [14°, 67°]. We don't show 
the region for 0.99 < /i < 1 because the plot descends steeply 
to zero in this region and it is difficult to show the nature of 
the plot of this region on this scale. 



We focused on the bending angle and the retarded time 
for a trajectory of the null ray with the impact parame- 
ter A which starts at some finite radius ro and propagates 
to the infinity. We obtained an analytical expression in 
terms of the elementary functions which uniformly ap- 
proximate the bending angle &\(ro) and the retarded 
time function T\(ro) with high accuracy. Knowledge of 
these two functions, $a(?"o) and T\(ro), allows one to cal- 
culate other quantities which are of interest in possible 
applications. 

For example, let us consider a scattering problem, 
when a photon with a given impact parameter A comes 
from infinity and after passing near the black hole goes 
to infinity again. The total bending angle for this process 
is 2$A(r*), where r* is the radius of the photon's turning 
point. This quantity can be approximated by 2$*, where 
is given by (146[) . Similarly, the expression ([71]) can be 
used for an approximation of the time delay quantities 



for the scattering problem. 

A knowledge of the bending angle function §\(r) is 
sufficient for a reconstruction a complete ray trajectory. 
Really, let us fix the impact parameter A and consider 
a ray with this impact parameter emitted at the radius 
ro in the D-dimensional Schwarzschild-Tangherlini space. 
Since the ray trajectories are planar, it is sufficient to con- 
sider rays lying in the equatorial plane passing through 
the point of emission. We choose cf>o = for this point. 
The ray trajectory is defined by the equation 

= F(r) = * A (r) - $ A (r ) • (92) 

To obtain an approximated form of the trajectory equa- 
tion it is sufficient to substitute <& a instead of the exact 
values 

Similarly one can use the approximating expression 
for solving the problem of the reconstruction of a null 
geodesic connecting two points (ro,(/>o) and (ri,0i). In 
the 4 dimensional case this problem can be reduced to 
one non-linear equation which contains as the arguments 
the Jacobi elliptic functions [35| . By using the developed 
approximation one can reduce this problem in any num- 
ber of dimensions to solving non-linear equations which 
contain only elementary functions. 

We hope that the obtained results will be useful in 
studying these and other problems connected with light 
propagation in the black hole vicinity which are of in- 
terest in the application to astrophysics. Since the devel- 
oped approximation is valid in a spacetime with arbitrary 
number of dimensions, it might be useful also for study 
black holes in the models with large extra dimensions. 
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